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We compare and contrast two types of deformations inspired by mixing applications - one from the mixing of fluids 
(stretching and folding), the other from the mixing of granular matter (cutting and shuffling). The connection between 
mechanics and dynamical systems is discussed in the context of the kinematics of deformation, emphasizing the equiv- 
alence between stretches and Lyapunov exponents. The stretching and folding motion exemplified by the baker's map 
is shown to give rise to a dynamical system with a positive Lyapunov exponent, the hallmark of chaotic mixing. On 
the other hand, cutting and shuffling does not stretch. When an interval exchange transformation is used as the basis 
for cutting and shuffling, we establish that all of the map's Lyapunov exponents are zero. Mixing, as quantified by the 
interfacial area per unit volume, is shown to be exponential when there is stretching and folding, but linear when there 
is only cutting and shuffling. We also discuss how a simple computational approach can discern stretching in discrete 
data. 



[. INTRODUCTION 

The essence of mixing of a fluid with itself can be under- 
stood in terms of an array of striations of, say, two different 
colors of the same fluid (or two different fluids such as cof- 
fee and cream) undergoing stretching and folding. On top of 
stretching and folding we may superimpose diffusion, reac- 
tion, and, in special circumstances, breakup processes lead- 
ing to droplet formation. 1,2 This approach is the backbone of 
lamellar models of mixing. A fundamental measure of the 
quality of mixing is ay, the interfacial area per unit volume of 
the striations (lamella or layers). Let S be the interfacial area 
between fluid layers within a volume V enclosing the point x 
at time f, then the interfacial area per unit volume is given by 2 

ay(x, f) = lim — . (1) 

v^o V 

A larger ay corresponds to better mixing. 

We can imagine many iterative mixing protocols that gen- 
erate large values of ay and create striations of the material of 
continually decreasing thickness in time. For fluids, a mul- 
titude of clever mixing designs can lead to the thinning of 
lamella, many inspired by a direct correspondence between 
the kinematics of mixing and chaotic dynamical systems. 3 
The simplest representation of mixing in terms of stretching 
and folding is the Smale horseshoe map, which stretches out 
a piece of material and folds it onto itself to form the shape 
of a horseshoe. A limiting case of this procedure is a map 
that stretches, cuts and re-stacks to generate interfacial area, 
the baker's transformation, named after the process by which 
a baker kneads dough. 

Granular mixing has been studied as well, but less exten- 
sively than fluid mixing. In many respects, the ideas applied 
to fluids carry over to granular matter. 4 A key difference be- 
tween the two is that granular flows may present surfaces of 
discontinuity, such as the interface between a flowing surface 



layer and the underlying static bed of granular material in an 
avalanche. 5 This new aspect of the flow leads to different mod- 
els for the kinematics. In particular, mixing in granular flows 
in rotating containers 6 ("tumblers") can be thought of as "cut- 
ting and shuffling," 7 a process different from stretching and 
folding. More on the fascinating behavior of granular matter 
can be found in Ref. |§l 

Stretching is a fundamental concept in mechanics and is 
covered in every continuum mechanics textbook in the con- 
text of kinematics (see, for example, the classic volumes of 
Truesdell 10 and Gurtin 11 ), where it is identified with shear or 
extensional strain. Cutting and shuffling, in contrast, has been 
explored only recently. To illustrate the fundamental differ- 
ence between the mixing mechanisms of stretching and fold- 
ing versus cutting and shuffling, we compare two types of sim- 
ple idealized mixing protocols: the well-known baker's map 
and cutting and shuffling maps based on interval exchange 
transformations. Along the way we discuss the elegant con- 
nection between continuum mechanics and dynamical sys- 
tems and show how some of the calculations of the relevant 
kinematic quantities, such as the deformation gradient and the 
principal stretches, are performed. Pertinent concepts from 
the theory of mixing are also reviewed within this context. 
We then discuss how to discern and measure stretching (or 
lack thereof) in practice. 



II. KINEMATICS OF DEFORMATION 

We restrict our discussiuon to two spatial dimensions for 
simplicity^ Consider a motion O from the undeformed (ref- 
erence or initial) configuration of the continuum (body) So to 
the deformed (current or final) configuration S, at time t , both 
of which are regions in the Euclidean plane. Keeping with the 
mathematics notation, in this paper, we denote sets by capi- 
tal calligraphic letters. Typically, we concern ourselves with 



2 



motions that map the body back into itself, that is, S t = So, 
but this restriction is not important for what follows. We are 
interested in discrete-time motions such as a repetitive mixing 
protocol. That is, we allow only t = nT, where T is the dura- 
tion (period) of the motion and n is a positive integer. Then, 
we may write the motion as a map 



<t>: 



(2) 



such that after one iteration every X in So is mapped to an 
x - O(X) in S t , where by x = (x, y) T we shall denote the po- 
sition vector in the deformed configuration and X = (X, Y) T 
represents the coordinates in the undeformed configuration. 
The T superscript denotes the transpose, meaning that (x,y) 
is a row vector and (x, y) T is a column vector. Unless other- 
wise noted, all vectors not written out in component form are 
considered to be column vectors. 

We can now define the deformation gradient (or Jacobian 
matrix of the map) as 



F = (Va;) T = 



d(x,y) 
d(X, Y) 



I dx dx\ 



By 



oy 



(3) 



ydX dY> 



It is typically assumed that the map is invertible and differ- 
entiable a sufficient number of times so that F exists and 
< det F < ooJJMi In this paper we relax the differentiability 
assumption to consider a wider (and, arguably, more interest- 
ing) class of motion. 

The polar decomposition theorem 2,10 ' 11 ' 13 allows us to write 



F = RU, 



(4) 



where U is a symmetric positive definite matrix 14 (due to 
the assumption that det F > 0) and R is a proper-orthogonal 
matrix. 15 That is, we can locally decompose the deformation 
into a rotation R and a stretch U. Because the matrix U is sym- 
metric positive definite, it has an orthonormal basis of eigen- 
vectors \e\,e2\ and strictly positive real eigenvalues {cr 1,0-2} 
satisfying 



Ue,- = <x,e,. 



(5) 



The eigenvalues cr, and eigenvectors e,- are called, respec- 
tively, the principal stretches and the principal directions be- 
cause an infinitesimal line segment of length dl oriented in the 
direction e, has length o-,d{ after undergoing the motion (that 
is, after the map O is applied). Hence, cr, > 1 for stretching 
and < cri < 1 for compression along the direction e,. 

Problem 1. Consider dx - 0(X + dX) - <S>{X), and show 
that dx = FdX+ higher order terms when \\dX\\ <sc 1. (Hint: 
use Taylor's theorem for a vector function.) Here, \\dX\\ = 
~\/dX ■ dX is the norm induced by the usual Euclidean dot 
(inner) product. Then let dX = (d£)ej and show that \\dX\\ = 
dl and \\dx^\ = cridl. 

The stretches may depend on the coordinate X in the un- 
deformed configuration and also on time (or, in the present 
context, the number of times n that the mixing protocol is re- 
peated). Typically, we are concerned with the largest principal 



stretch & = max{cri, aj\, which is a scalar field & = &{X; n) 
that describes the stretching experienced by the body due to its 
motion. In practice, it is convenient to calculate the eigenval- 
ues {ki , K2} of the (right) Cauchy-Green strain tensor C = F T F 
instead. These are just the squares of the principal stretches. 
In this way & can be computed without explicitly finding the 
polar decomposition of F. 

Problem 2. Show that k,- = cr 2 using the definition of C, the 
polar decomposition theorem, and the properties of eigenval- 
ues. 



III. DYNAMICAL SYSTEMS FRAMEWORK OF KINEMATICS 

Equation (f2]i also defines a dynamical system, which in 
the most general sense is defined as a rule of evolution on 
a state space (the body). 16 This connection between the kine- 
matics of continua and dynamical systems has been success- 
fully exploited in the study of both fluid mixing 2 and granular 
mixingii 

There is a direct correspondence between the languages of 
dynamical systems and continuum mechanics, the most im- 
portant of which, for the present purposes, is the correspon- 
dence between stretches and Lyapunov (characteristic) expo- 
nents. For a discrete-time map, these exponents are defined 
(see, for example, Ref. \ J 3 Sec. 5.3.1) as 



A(X,v) = lim -ln||(VO") T u|| 



(6) 



where <[>" = O o • • ■ o <S> (n compositions of the map). The 
Lyapunov exponents depend on the position X and on the 
direction v Q\v\\ = 1). The quantity VO" can be calculated 
by the chain rule along a trajectory starting at a given X in 
So- If F = (VO) T happens to be independent of X, we have 
(VO") T = F". Lyapunov exponents are important in the con- 
text of the asymptotic stability of infinitesimal perturbations 
and, provided certain conditions are satisfied, can be inter- 
preted as the growth (or decay) rates of these perturbations 
along a trajectory. 

A related concept is fixe, finite-time Lyapunov exponents de- 
fined as 



y(X,v;n) = - In ||(VO") T v|| 
n 

Note that if (VO»") T = F", we have 



(7) 



max 

v±0, ||w|[=l 



y(X, v,n)=^ In ^p((F") T F"), (8) 



where p(A) denotes the spectral radius of the matrix A, that 
is, its largest eigenvalue in absolute value. As we discussed 
at the end of Sec.|lll the square root of the largest eigenvalue 
of (F") T F" is the largest principal stretch (for the deformation 
resulting from applying the map n times). Therefore, there is a 
one-to-one relation between the largest finite-time Lyapunov 
exponent and the largest principal stretch: 



y max (X; n) = max y(X, v; ri) 
v*0, |M|=1 



1 



lno-(X;«). (9) 
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If the limit in Eq. © exists, then this relation carries over to 
the largest (infinite-time) Lyapunov exponent as well. 18 

Problem 3. Derive Eq. (0 using the Rayleigh-Ritz theorem 
from linear algebra. 1 - Hint: write out the norm in the defini- 
tion of y as a dot product and re-arrange terms. Also, you are 
allowed to bring the max inside any logarithm or square root 
because these functions are monotonically increasing. 

We will call a dynamical system with a positive Lyapunov 
exponent chaotic (although there are variety of sometimes 
equivalent ways to define chaos 1 -). If a motion in the sense 
of continuum mechanics stretches the continuum (<x > 1), 
the corresponding dynamical system has a positive Lyapunov 
exponent by virtue of Eq. ((9). Note that, in this context, 
the chaotic dynamics occur in physical space, not in phase 
space as is the case for the dynamical systems of classical 
mechanics^ 

There is a temptation to equate stretching and folding mo- 
tions (equivalently, ones that give rise to chaotic dynamical 
systems) with efficient generation of ay. However, as we will 
show in the following examples, there are other possibilities 
for generating ay without stretching and folding. 



cally the resulting rectangle in half along X - 1/2 and stack- 
ing the pieces [see Fig. Q] (BM)] - much like how a baker 
kneads dough or a taffy pull machine makes candy. 26 Other 
applications include the Kenics mixer— and the related par- 
titioned pipe mixer (see, for example, Ref. 0, Sec. 8.2). The 
baker's map has even been used to explain the movement of 
bubbles in a foam network. 28 In fluid mixing, the cut and re- 
stack step cannot be accomplished exactly, but cutting and re- 
stacking is precisely what happens in the extrusion of multi- 
layer plastics^ 

Even though the baker's map constitutes stretching, cutting 
and re-stacking, it possesses (as we will show) the essential 
property of a stretching and folding motion - stretching in 
the continuum mechanics sense (& > 1) or, equivalently, a 
positive Lyapunov exponent. The Smale horseshoe map 2,30 ' 31 
is classical example of stretching and folding. However, the 
horseshoe map does not preserve area, which has significant 
implications for the types of chaos exhibited by the dynamical 
system because it allows for the existence of attractorsM We 
avoid such complications by idealizing the Smale horseshoe 
map as the baker's map. 

From Eqs. <[3j and ( fTOb it follows that 



IV. EXAMPLE MIXING PROTOCOLS 

One of the best known stretching and folding operations is 
the baker's map, 2,17 which can be found as early as 195 1 in the 
mixing literature 21 and 1937 in the mathematics literature. 22 
In contrast, cutting and shuffling, two examples of which we 
construct using interval exchange transformations, is a type of 
mixing protocol that has been proposed only recentlyi^ 2 ^— 
Suppose we begin with a square region, half of which is filled 
with a while material while the other half is filled with a gray 
material (see Fig. [TJ. Then, from an inspection of the initial 
(undeformed) and final (deformed) configurations of the con- 
tinuum undergoing deformation described by these protocols, 
it is clear that two of them — the baker's map (BM) and the 
second cutting and shuffling map (CS2) — lead to mixing of 
the gray and white materials. (The map CS1 does not pro- 
duce further striations of material, as it only re-arranges the 
layers produced by the first set of cuts in a periodic manner.) 
However, the following simple analytic calculations, which 
are confirmed by numerical calculations in Sec.[Vl] show that 
only the baker's map stretches the material continuum and 
gives rise to a chaotic dynamical system in the usual sense. 



A. The baker's map 

We let S, = Sq = [0, l] 2 , and write the baker's map B m: 
[0, l] 2 -> [0, l] 2 as 



O bm (X) 



2X,-Y , 0<X<- 

i IV i (10) 



Simply put, this map involves compressing the unit square to 
half its height, stretching it to twice its width, cutting verti- 



Frm — 



2 



X + 



1 



(11) 



and Fbm is undefined at the cut X = 1/2. The polar decom- 
position of this matrix is simple: R BM = I and Ubm = Fbm- 
The eigenvalues of Ubm are en = 2 and o~2 = 1/2. Therefore, 
the largest principal stretch is & — 2 > 1 everywhere after 
one iteration of the protocol (n = 1). Clearly, the baker's map 
stretches the underlying material continuum. 

Because F BM is independent of X, iterating the baker's 
map results in 



(v*bm) T 



BM 



2 



2" 
2" J 



(12) 

and Fbm is undefined along the cuts, that is, the set of points 
C" M = LY in [0,l] 2 |X = //2", i=l,2, 



Once again, R B 



BM 



I and U" BM 



BM' 



.,2"-l; < Y < 1} 
(13) 

Consequently, the 



largest principal stretch of the deformation is 
6- BM (X;n) = 2 n (X £ C BM ) 



(14) 



Turning to the Lyapunov exponents, the two directions of 
interest are the principal directions e, (the eigenvectors of 
Ubm), which are e\ = (1,0) T and e 2 = (0, 1) T . Therefore, 
upon multiplying the vectors e\, e 2 by the matrix in Eq. ( fT2b 
and taking the norm of the resulting vector, In ||Fg M ei,2l| = 
+n In 2. We substitute this expression into Eq. (0 and obtain 
the finite-time Lyapunov exponents (along the principal direc- 
tions) of the baker's map 

y BM (X, ei , 2 ;») = ±ln2 (X^Q, (15) 

and they do not exist for X in Cg M . 
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Problem 4. Show that the larger finite-time Lyapunov ex- 
ponent in Eq. ( fT3T > can also be calculated using Eqs. (0 and 

CD. 

The Lyapunov exponents defined in Eq. © can be obtained 
by taking the limit n — > oo in Eq. ( fT~5T > : 

A B M(X > ei )2 )=±lii2 (X^Q. (16) 

A subtle but important point 33 is that the set of cuts Cg M is a 
null subset of the domain So as « — > oo 31 because it is a set of 
measure zero. Thus the quantities in Eqs. (TT~5T > and ( TTol l exist 
almost everywhere. 

Since /l BM (X, ei) > 0, it follows that the baker's map gives 
rise to a chaotic dynamical system; that is, if we were to track 
an infinitesimal ball of material points of the continuum, they 
will spread exponentially fast from each other with repeated 
applications of the map because it possesses a positive Lya- 
punov exponent. The simplicity of the baker's map and its 
ability to render many analytical calculations tractable have 
made it one of the classic examples of a chaotic dynamical 
system. 34 We refer the reader to, for example, Ref. 3_5J Sec. 5.1 
for a thorough overview, including the non-area-preserving 
(dissipative) version of the map. 

B. Cutting and shuffling maps 

Cutting and shuffling is an operation well-known to card 
players and mathematicians. One of the simplest dynamical 
systems that can be successfully studied analytically is the 
interval exchange transformation; see, for example, Ref. 
Sec. 14.5. An interval exchange transformation, which we 
write as Esji'. I — * I, replaces a deck of cards with a con- 
tinuous interval I of the real line (we take I = [0, 1] with- 
out loss of generality) and subdivides it into a collection of 
k disjoint subintervals S = . . . ,Ik}- These subintervals 
are translated ("shuffled") according to a rule given by some 
permutation II of the integers between 1 and k. Finally, the 
interval I is put back together as I — I~na) U ■ ■ • U In(t)- The 
unusual type of continuum motions encountered in tumbled 
granular flows^ 2 ^— that we noted in Sec. I are specific two- 
dimensional generalizations, called piecewise isometries } 1 ^ 
of these simple maps. 

To mix the domain £>o = [0, l] 2 , we can apply an inter- 
val exchange transformation in the T-direction and extend it 
in the X-direction by making each subinterval I~, into a rect- 
angle of unit horizontal length. This construction is a spe- 
cial case of the more general class of rectangle exchange 
transformations} 9 Consider the simple special case of k — 
8 equal subintervals shuffled according to the permutation 
II([12345678]) = [15263748], numbering the subintervals 
from bottom to top so that intervals 1 through 4 are gray and 
intervals 5 through 8 are white. This interval exchange trans- 
formation is equivalent to cutting a deck of 8 cards exactly 
in half and shuffling them perfectly so that the bottom card 
from the first half (gray) is just below the bottom hard from 
the second half (white), and so on. As shown in the middle 
path of Fig. [T](CS1), one iteration of this map results in lay- 
ers of gray and white that appear identical to those produced 



by two iterations of the baker's map. It is not difficult to see 
how one iteration of the same cutting and shuffling map with, 
say, k — 64 subintervals results in the same number and place- 
ment of layers of gray and white as n — 5 iterations of the 
baker's map. However, repeated application of the map CS 1 
does not result in mixing for any choice of k. At n — 2 (for 
k = 8), demixing 40 occurs: the number of distinct stripes de- 
creases from 8 to 4. Fortunately, there is a well-developed 
theory to guide us in constructing interval exchange transfor- 
mations that mix well. 

To this end, consider the more general case in which the 
intervals are not equal and are not shuffled in such an or- 
derly fashion. It is well known that for an interval exchange 
transformation to exhibit interesting behavior (specifically, to 
have no periodic orbits and to be ergodic 41 ), it must satisfy the 
Keane condition. 42 This condition requires that the ratio of the 
lengths of adjacent intervals |J,-|/[Ji_i| be an irrational number 
and that the permutation II be irreducible; that is, applying II 
to any of the subsets {1}, {1, 2), {1, 2, 3} up to {1, 2, . . . , k - 1} 
does not yield a permutation of just the elements of the sub- 
set. For example, the permutation II([12345]) = [31254] is 
reducible (that is, not irreducible) because the first three ele- 
ments are a permutation of only themselves (neither 1 , 2, nor 
3 maps to 4 or 5). In contrast, II([12345]) = [31524] is irre- 
ducible because 3 maps to 5. To satisfy the Keane condition, 
we consider k — 4 subintervals, and suppose their lengths 
are chosen so that the first is - rj and each consecutive 
subinterval has length r times the length of the previous one: 
= r|_T;_i|. Because the length of I — [0, 1] must be pre- 
served, rj + rrj + r 2 rj + Ptj — I. Given r we can solve for rj from 
this relation. We take r = VO and, numbering from bottom 
to top, Il([1234]) = [2431] to satisfy the Keane condition. 
In terms of the continuum mechanics language introduced in 
Sec. |nl the map takes the form 

®CS2(X) = 

\X,Y + (r + r 2 + P)j]) T , (0<Y<T]) 
(X,Y-tj) t , (77 < Y < [1 + r]rf) 

' (X, Y + (r 3 - 1)?7) T , ([1 + r]j] < Y < [1 + r + r 2 ]i]) 

(X, Y - (1 + r 2 )rj) T , ([1 + r + r 2 \r] < Y < 1). 

(17) 

The action of this cutting and shuffling map, constructed from 
an interval exchange transformation, is illustrated by the right 
pathinFig.[T](CS2). 

Problem 5. Derive Eq. dT7T > by considering where the rect- 
angles go in Fig.[T](CS2). (Hint: rectangle 1 moves up by the 
height of rectangles 2, 3, and 4, that is, (r + r 2 + r 3 )rj, etc.) 

Because we are simply translating the strips in the vertical 
direction, calculating the deformation gradient from Eq. (01 is 
trivial: 

Fcs2 = (J J) (XtC l CS2 ), (18) 

and it is undefined along the cuts, that is, the set of points 
C l CS2 = {X in [0, 1] 2 |0 < X < 1; Y = f-q, i = 0,1,2). 
Unlike for the baker's map, finding the set of cuts C" C51 for any 
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FIG. 1. Illustration of mixing by three different protocols. The left path (BM) shows two applications of the baker's map, which stretches the 
underlying continuum. The middle path (CS1) is one iteration of a cutting and shuffling map with k = 8 intervals and length ratio r = 1, which 
does not stretch and leads to the same macroscopic picture as the left path (BM). The right path (CS2) is two applications of a cutting shuffling 
map based on an interval exchange transformation with k = 4 and r = VL3, which also does not stretch but results in a different macroscopic 
picture. Only the maps BM and CS2 produce an increasing number of finer striations of each color upon repeated application. 



n for the above cutting and shuffling map (CS2) is not trivial 
and, in general, an open problem. 43 Still, it is a countable set, 
and therefore C^ S2 is a null subset of Sq. 

Problem 6. Write a computer program that finds the number 
and location of the distinct cuts after n applications of Ocs2- 
What happens if you change r or FI? 

For any given trajectory that does not start at some X in 
Cq S2 , it is clear that the deformation gradient after n iterations 
is the product of identity matrices like the one in Eq. (flSt , 
Hence F£ S9 = I. Therefore, the principal stretches are o~\ = 
0"2 = 1, and 

<T CS 2(X;«)=1 (XtC^). (19) 

From Eq. (O and F£ S2 = I, we can calculate the finite-time 
Lyapunov exponents along the principal directions: 

y CS2 (X,e 1 ;n) = y CS2 (X,e 2 ;n) = (X<tC n CS2 ), (20) 

and they do not exist along the cuts. If we take the limit n — > 
oo in Eq. d20l >. we arrive at the Lyapunov exponents [Eq. ©]: 

Ac S2 (X,e l ) = A CS2 (X,e 2 ) = (X * C^)- (21) 

Equations d20"l > and (|2"TT > are in stark contrast with Eqs. (Q3) 
and ( fT6b . In addition, this type of argument can be formalized 



into a proof showing that all Lyapunov exponents are zero 
for the more general class of piecewise isometry maps on the 
plane— 

Problem 7. Show that the finite-time Lyapunov exponents 
in Eq. d20b can also be obtained from Eqs. © and (fT9l . 

Therefore, the map Ocs2 (and, by the same logic, the map 
<I>csi) does not give rise to chaotic dynamics. Specifically, 
an infinitesimal ball of material points in the continuum does 
not spread apart under n iterations of this map (for any choice 
of «), unless the ball happens to overlay a cut. Hence, the 
distance between any two points in So can grow only if they 
become separated by a cut. 

We have assumed that the same specific shuffle (defined by 
the given permutation n) is performed once at each iteration 
of the map. In doing so, we did not consider the possibility 
that a white piece of material may be moved next to another 
white one, meaning that the shuffle is not necessarily optimal 
in terms of mixing. Thus, an improvement would be to allow 
multiple shuffles to be performed at each step of the proto- 
col. There is an elegant theory 45,46 (in the context of shuffling 
decks of cards) that gives an estimate of how many such shuf- 
fles are required for the stack to be "sufficiently random." 
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V. DEFINING AND QUANTIFYING MIXING 

The goal of the protocols we have considered is to mix the 
gray with the white in Fig.Q] but what do we mean by "mix"? 
If we go by the earlier intuitive definition, that is, to mix is to 
generate interfacial area ay, then one thing is clear about the 
baker's map protocol from Sec. |IV A| the number of stripes 
(of unit length) in the final configuration is 2" +1 . Thus, the 
bulk interfacial area between the gray and white materials is 

fl V = 2" +1 -l =2e ta -l. (22) 

By expressing a y as an exponential function of n with growth 
rate h = In 2 = yj^, 47 we are able to illustrate the well-known 
result that chaotic advection, through stretching and folding, 
generates ay at an exponential rate. 

For the cutting and shuffling maps from Sec. lIVBl the most 
interfacial area that can be generated in one iteration is equal 
to the number of cuts, which is k — 1 for k subintervals. Hence, 
assuming that the interface at Y = 1 /2 in the initial configura- 
tion does not coincide with a cut, the upper bound is 



a v < 1 + (k - l)n. 



(23) 



Equation ( l23l is only an upper bound because not every cut 
gives rise to an intermaterial interface after shuffling. For ex- 
ample, while on the first iteration of the map CS1 in Fig.[T]the 
seven cuts produce seven interfaces, on the second iteration 
only three interfaces remain though we have made seven fresh 
cuts. The nature of the upper bound in Eq. d23l illustrates 
a general principle about interval exchange transformations 
and a conjecture about their cousins piecewise isometries: the 
growth of ay is sub-exponential (specifically, algebraic) be- 
cause of the lack of stretching by the map. 

Another way to define mixing is by using ideas from the 
ergodic theory of dynamical systems. There, (strong) mixing 
is succinctly defined (see, for example, Ref. 17, Sec. 3.7 or 
Ref.® as 

lim Area(O n (^i) n^b) = Area(yb) Area(^ 2 ), (24) 

n— >oo 

where yii and ^2 are any two subsets of So and n denotes the 
intersection of two sets (that is, the material they have in com- 
mon). Put simply, under iteration of the map O, the set J{\ 
eventually becomes spread out evenly throughout the domain 
So so that no matter what other set ^2 we choose, the amount 
of material in it that came from J[\ is the same. In a more for- 
mal definition, Area( ) is replaced by the appropriate invariant 
measure. The maps we consider here are area-preserving, and 
therefore area is the proper measure of "size" for the defini- 
tion. 

Following the argument of Ref. [l_7| (Theorem 3.7.2), sup- 
pose J{\ is the gray region of the initial configuration So (re- 
call Fig.[T}, that is, S{\ = [0, 1] X [o, 5]. Then, no matter what 
set we pick for ^2, the number of stripes of each material 
within this volume grows exponentially, eventually resulting 
in half of each. Therefore, as n — > 00, Area (<&"(j?li) n M2) — > 
\ Area(yt 2 ). Because Area(yti) = \, it follows that the left- 



and right-hand sides in Eq. ( l24b are equal, and the baker's map 
represents strong mixing. 

For the cutting and shuffling maps it should be clear that 
this argument breaks down. The second cutting and shuffling 
map (CS2) in Fig.Q]illustrates how the gray and white stripes 
under this map have non-uniform thickness and a pair of pre- 
viously cut gray pieces may be "glued back together" at a later 
iteration. Thus, a cutting and shuffling map based on an inter- 
val exchange transformation can only be shown to be weakly 
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mixing 

1 n_1 

lim - V I AreatO'^O n Jl 2 ) - Area^O Area(^ 2 )| = 0. 

1=0 

(25) 

This proof is far too involved to present here. Put simply, 
Eq. d25l > relaxes the requirement that the same amount of ma- 
terial from the subset J[\ be found in any other subset ^2. 
Instead, it is only required that this be true on average, over 
many iterations. Still, weak mixing is a stronger result than 
the ergodicity property mentioned previously that an inter- 
val exchange transformation acquires by satisfying the Keane 
condition^ 

To summarize, according to both the practical measure of 
mixing defined in Eq. ([T]) and the mathematical measures of 
mixing defined in Eqs. (124-b and (1251 . both of these protocols 
mix, albeit at different rates and with different "strengths." In 
this respect, the cutting and shuffling map defies conventional 
wisdom - there is no stretching, all Lyapunov exponents are 
zero, and the trajectories of material points in the continuum 
under this map are not chaotic. 51 



VI. INFERRING STRETCHING FROM DISCRETE DATA 

Suppose we performed an experiment with white and gray 
putty and recorded the initial and final configurations. We 
would like to find out how the putty was stretched due to the 
mixing protocols depicted in Fig. Q] In a typical experiment 
only a finite number of material points in a continuum can be 
tracked. In other words, some collection of M 2 material points 
{Xjj}i j=\ M in So in the undeformed configuration are iden- 
tified, and their locations {a; ; ,j},, ; =i...M in S t=n T in various de- 
formed configurations (for example, for different n) recorded, 
as shown schematically in Fig. [2] From these values the mo- 
tion and deformation of the continuum can be reconstructed. 

A particularly simple and effective numerical approach to 
achieving this reconstruction is the standard central difference 
approximation to the deformation gradient [recall Eq. fl3)]: 



F(X y ;n) = 
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This approximation is used, for example, in extracting La- 
grangian coherent structures (barriers to mixing and trans- 
port) from experimental and simulated data, 53 in determining 
the type (hyperbolic, elliptic, or parabolic) of periodic points 
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(that is, X such that <i> p {X) = X, p = 1,2, . . .) of mixing 
protocolsri^ and in calculating the largest principal stretch 
in fluid mixing experiments with laminar 56 - 57 or turbulen t 58 ' 59 
velocity fields. Understanding stretch fields is of immense 
practical importance in both the industrial and laboratory set- 
ting from the micro to the planetary scale. Their use in finding 
barriers to transport and mixing in the ocean was discussed in 
RefM 




FIG. 2. Cumulative movement of material points of the undeformed 
continuum S after n iterations of O, resulting in the deformed con- 
figuration B,=„T- 

Tracking the motion of a collection of material points, as 
illustrated schematically in Fig. [2] is the objective of this nu- 
merical approach. It may seem that many rearrangements of 
the continuum can result in a non-trivial deformation gradient, 
at least somewhere in the domain. Therefore, it would appear 
that stretching occurs under most maps. However, such intu- 
ition can often be wrong. 

Figure |3] shows the numerical results for & based on track- 
ing the movement of M 2 = 251 2 uniformly-distributed mate- 
rial points under each of the maps in Fig. Q]for the number 
of iterations of each map depicted there. Note that the three 
narrow vertical lines of high & values for the baker's map and 
similar horizontal lines for the cutting and shuffling maps are 
due to the cuts in these protocols causing the stretch to be in- 
finite there (that is, F" is undefined). The lines are vertical 
for the baker's map, because the cuts are vertical. The middle 
line in Fig. [2 a) corresponds to the cut that is made in the first 
application of the map, and the other two lines correspond to 
the cut made in the second application of the map, which il- 
lustrates the set Cg M defined in Eq. ( fT3l ). For the cutting and 
shuffling maps, the lines are horizontal because the cuts are 
horizontal. Specifically, the dark horizontal lines in Fig. [2b) 
depicts the set of cuts C X CSI = {X in [0, l] 2 1 < X < 1; Y = 
//8, i = 1, . . . , 7), while those in Fig. [2c) depict the set of 
cuts C£ S2 (n = 2) noted after Eq. ( fT~8b . 

The apparent thickness of these lines is due to the coarse- 
grained view of the deformation that we obtain from the nu- 
merical data. Taking more points (larger M) makes these lines 



arbitrarily thin. (Their thickness is w 2/M, that is, twice the 
spacing between points because we used a central difference 
approximation to F.) If we had applied the map more times 
(larger n) and/or had fewer material points, then the resulting 
thicker cuts would completely obscure the picture. 

In Fig. [2 the expected results are found. We applied the 
baker's map twice [recall Fig.Q](BM)], so that n — 2. From 
Eq. ( TBI , the largest principal stretch is & — 2 2 = 4. We can 
also calculate & from the numerical data by taking the square 
root of the largest eigenvalue of [Fbm(^; 2)] t Fbm(^; 2) (re- 
call the discussion at the end of Sec. [TTJ, where F is defined 
in Eq. (|26T i. In complete agreement with the theoretical result, 
the numerical calculation produces a "background" value (that 
is, the value of & away from the cuts discussed above) of 4, 
as shown in Fig. [2 a )- This result can also be anticipated by 
realizing that the baker's map stretches the continuum in the 
X-direction by a factor of 2. Applying the map twice gives a 
stretch ratio of 4. Similarly, for both of the cutting and shuf- 
fling maps, the largest principal stretch is equal to 1 for all n 
from Eq. ( fT9] l, that is, the maps do not stretch. The numerical 
calculation shown in Figs. Ob) and (c) confirms that & = 1, 
In both cases, the background value in the plots is 1, and & 
is different from the "background" value only along the cuts, 
where it is infinite. 



VII. CONCLUSION 

The mixing of continua can be accomplished by a great 
variety of maps of different complexities. In mechanics, the 
central theme is stretching. Stretching leads to a positive Lya- 
punov exponent and chaos when the problem is translated into 
the language of dynamical systems. For mixing, the central 
theme is generation of interfacial area per unit volume ay. 
Stretching and folding (equivalently, chaotic dynamics) is an 
efficient way to do this. 

Surprisingly, cutting and shuffling can also "rearrange" ma- 
terial points in a continuum quite well. The cutting and shuf- 
fling maps, which would be considered pathological in "clas- 
sical" continuum mechanicsrSiii do not stretch, possess no 
positive Lyapunov exponents, and exhibit no chaotic behavior 
in the usual sense, yet they mix. Even though interfacial area 
is produced at an asymptotically slower rate than for a chaotic 
motion that stretches, for short times it appears possible for 
cutting and shuffling to dominate stretching and folding, de- 
pending on the number of cuts k in Eq. ( l23l l and the growth 
rate h in Eq. d22b . Although we may be misled to conclude 
that a cutting and shuffling map stretches the underlying con- 
tinuum based on the complicated pattern of non-uniform stri- 
ations produced by it, analytical and sufficiently-refined nu- 
merical calculations show otherwise. 

Recent work- 7 ^ 2 ^ 5 . takes the idea of cutting and shuffling 
even further, arguing that it provides the "skeleton" of certain 
regimes of granular flow in tumblers. In the present work, 
we illustrated how cutting and shuffling can also lead to mix- 
ing without chaos. And, in more complicated maps, cutting 
and shuffling can lead to complex dynamics without chaos in 
any usual sense of the wordier— Thus, cutting and shuffling 
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(a) left path in Fig. 1 (BM) (b) middle path in Fig. 1 (CS1) (c) right path in Fig. 1 (CS2) 



FIG. 3. Largest numerically-computed principal stretch field &(X\n) for the protocols depicted in Fig.Q] The background value is 4 in (a) but 
equals 1 in (b) and (c). 



maps and their generalizations are more than just a mathemat- 
ical exercise, and have the potential to accurately describe the 
underlying framework of mixing in physical systems. 
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